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Abstract 

Why sex evolved and it prevails in nature remains one of the great 
puzzles of evolution. Most biologists would explain that it promotes ge- 
netic variability, however this explanation suffers from several difficulties. 
What advantages might sex confer? The present communication aims at 
certain investigations related to this question, in this way we introduce 
sexual recombination on the Hotzel model (with males and females) and 
we compare these results with those from asexual reproduction without 
recombination. 
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1 Introduction 

Sex, which involves the alternation of meiosis and gamete fusion, is a rather in- 
efficient means of self-propagation as compared to asexual reproduction, where 
offspring stem only from a mitotically produced cells. One of the most com- 
mon reasons used to explain the origin and maintenance of sex is its ability to 
reduce the mutation load if consecutive mutations lead to an increasing decline 
in relative fitness, although it is not clear a priori that the heritable variance in 
fitness is significantly increased by sex. Despite decades of developing theoreti- 
cal models to explain why sex is such a widespread phenomenon and how sexual 
reproduction may confer advantages that outweigh its disadvantages, until now 
no such general clear advantage has been found. 

Investigations of evolutionary problems by physicists have in fact boomed in 
the last few years. Since computer simulations of natural systems can provide 
much insight into their fundamental mechanisms, they can be used to test the- 
oretical ideas that could be otherwise viewed as too vague to deserve the status 
of scientific knowledge [1, 2]. In this way, many computer models in population 
dynamics have been proposed to investigate the evolution of sex and its justi- 
fication, as well as the comparison between sexual and asexual reproduction, 
for instance, the Redfield model [3], the Penna bit-string model [4], a genomic 
bit-string model without aging [5] and Stauffer model [6, 7, 8]. 
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Of particular interest here is the Heumann-Hotzel model [9] , which originally 
simulated the evolution of asexual population, composed of haploid individuals, 
without recombination. Thus now wc introduce the recombination in this model, 
in order to find out if the sexual reproduction (with males and females) can 
produce better results than the simple asexual reproduction [10]. In the next 
section, we describe the standard and the modified Hcumann-Hotzel model, in 
section 3, we present our results and in section 4, our conclusions. 

2 Heumann-Hotzel Model 

Since it has been proposed by Michael Heumann and Michael Hotzel in 1995, 
the Heumann-Hotzel model [9], which was an unsuccessful attempt to intro- 
duce more ages in the Dasgupta model [11], has remained forgotten due to the 
fact that after many generations it reduces to the two-age model of Partridge- 
Barton [12]. The Dasgupta model consists in taking into account some modifica- 
tions such as hereditary mutations and food restrictions in the Partridge-Barton 
model. In fact, the Heumann-Hotzel paper [9], treats basically the computer 
simulations using the generalized Dasgupta aging model proposed by Michael 
Hotzel in his work, imdcr the supervision of Dietrich Stauffer, in order to obtain 
the license to teach in German secondary school [13]. Michael Heumann, who 
was another teacher's candidate, worked only on the inclusion of the "Dauer" 
state in the Dasgupta model [14]. 

Recently, the Heumann-Hotzel model was reinvestigated and, according to 
the authors, with "simple and minor change in the original model" this inca- 
pacity to describe populations with many ages seems to be surmounted [15]. 

In the original version of the Heumann-Hotzel, the genome of each (hap- 
loid) individual is represented by one set of probabilities po; Pi )P2, •••,Pmaxage-i; 
where pa is the survival probability that an individual has to reach age a + 1 
from age a. At every time step t, iV(t) *maxage individuals are chosen randomly 
to have their survival probability Pa altered by mutations to p'a = Pa* exp(e), 
where the age a is also randomly chosen. N{t) is the size of the population at 
time t and maxage is the maximum age one individual can live, which is set up 
in the beginning of the simulation. The quantity e is chosen randomly as any 
number between ei and £2 and when it is negative (positive) it corresponds to 
a deleterious (beneficial) mutation. 

The effect of food and space restrictions is taken account by an age- indepen- 
dent Verhulst factor, which gives to each individual a probability [1 — A^(t)/A^max] 
of staying alive; N^^^^ represents the maximum possible size of the population. 
This mean-field probability of death for the computer simulations has the benefit 
of limiting the size of population to be dealt with. The passage of time is 
represented by the reading of a new locus in the genome of each individual 
in the population, and the increase of its age by 1. After taking account the 
natural selection and the action of Verhulst dagger, at the completion of each 
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period of life, each individual gives birth to one baby (age=0) which inherits its 
set of probabilities (po,P'nP2i -vPrnaxage-i)- 

In the recent reinvestigation of this modeI[15], individuals with age a in the 
interval amin < a < ttmax will generate h offspring and the mutations are allowed 
only on a fraction F(0<i^<l)of the babies. 
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Figure 1: Total population for asexual reproduction (squares), sexual reproduc- 
tion: case (a) (triangles) and case (b) (stars). The inset shows the corresponding 
survival probability. 

In the sexual version, each (diploid) individual of the population, which 

consists of males and females, is genetically represented now by two sets of 
survival probabilities, -P(a)^ and -P(a)^, to be read in parallel. In this way, we 
have studied the following cases (see below): 



P{af = {vl,Vi,v\, ■■■,Vn 



maxagc— 1/ 



P{af = (Po>Pl>P2,->Pmaxage-l) 
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• Case (a) - The effective survival probability in some age will be the arithmetic 
average of the values present in both sets at that age: 

P/'/1^«ffective _ /Po+Po Pi + Pi Pmaxage-1 + Pmaxage-1 

na) 2 ' 2 2 ' 

• Case (b) - The effective survival probability in some age will be the maximum 
value between the values present in both sets at that age: 




Figure 2: Survival probability for asexual reproduction (circles), sexual repro- 
duction with ttmax = 17 and amin = 1 (triangles), 3 (stars), 5 (squares). 

If the female succeeds in surviving until the minimum reproduction age amin, 
it chooses, at random, an able male to mate (amin < age < Omax) and it gener- 
ates, with probability pi, b offspring every iteration until the maximum age of 
reproduction Omax- The offspring inherits its set of survival probabilities from its 
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Figure 3: Survival probability for asexual reproduction (circles), sexual repro- 
duction with flmax = 17 and amin = 1 (triangles), 3 (stars), 5 (squares). 

parents in the following way: the two sets of survival probabilities of the male, 
for instance, arc broken in the same random position, and the complementary 
pieces, belonging to different strings, are joined to form two male gametes. One 
of the gametes is then randomly chosen to be passed to the offspring. After that, 
rrim random mutations are introduced into this gamete, and the final result cor- 
responds to one string of the baby genome. The same process occurs with the 
female genome, generating the second string of the baby, with m/ mutations. 
At the end the offspring genome contains a total of M = mm + mf mutations. 
Finally, the sex of the baby is randomly chosen, each one with probability 50%. 
This procedure is repeated for each of the b offspring. 
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3 Results 



The simulation starts with individuals (half for each sex) and runs for 
400, 000 time steps, at the end of which (the last 10, 000 steps, when the popu- 
lation was stabilized) averages are taken over the population. The parameters 
of the simulations are: 

• Initial population A^o = 10, 000; 

• Maximum population size A'max = 100 * No; 

• Probability to give birth pb = 1.0; 

• Birth rate b = 1.0; 

• Mutation rate nim = ™/ = 1 per gamete; 

• ei = 0.02 and €2 = —0.04 (the same values used in the original Hotzel model). 

Our figures with = 10, 000 are confirmed by larger simulations with A^o = 
100, 000, and also by larger simulations with lO'^ time steps. 

From Figure 1 and its inset we can see that the diploid sexual population is 
not only larger than the haploid asexual one, but also presents a higher survival 
probability. 

In Figure 2 (case (a)) and Figure 3 (case (b)), we present the survival proba- 
bility as a function of age for different period of reproduction (amin < i < flmax)- 
Aging starts with reproduction: the survival rate decays as soon as reproduction 
age is reached. There are no individuals alive older than the maximum repro- 
duction age Cmax- Figure 4 corresponds to the case in which all the individuals 
of the population reproduces only once - the so-called catastrophic senescence 
effect [12, 16]. In this way, two rules of reproducing were adopted: 1) The 
reproduction age is the same for all individuals (omax = ctmin = 10), 2) The 
reproduction age is randomly chosen between amin = 5 and Omax = 10. We can 
noticed that this effect is more pronounced for the former [12, 16], as well that 
the responsible for that are both breeding once and breeding for all individuals 
at the same age. The explanation for these effects observed in Figures 2-4 is 
based on the Darwin theory: individuals must stay alive in order to reproduce 
and perpetuate the species. If they can no longer generate offspring but remain 
in the population, they are killed by the accumulation of deleterious mutations 

[!]■ 

In fact, real mutations can be divided into the common recessive (almost 
90% of the real mutations are recessive) and the rare dominant mutations. In 
this way, if among the many genes of a species, one of the father's genes differs 
from the corresponding one of the mother, then it adversely affects the child 
only if the mutation is dominant. Recessive mutations affect the child only 
if both the father and mother have them. In order to take into account this 
aspect in the sexual version of Hotzel, at the beginning of the simulation we 
choose randomly d dominant positions and keep them fixed during the whole 
simulation. The effective survival probability in the dominant positions (ages) 
will be the smallest value of the two located in the same position in both strands, 
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Figure 4: Survival probability for sexual reproduction in case of reproduction 
only at the reproduction age Omax = dmin = 10 (circles), and in case of repro- 
duction at some age between amin = 5 and amax = 10 (triangles). 

and for recessives ones the effective survival probability will be the arithmetic 
average of them. In Figure 5, we can see that the inclusion of dominance 
does not alter the lifespan of the population, although it has been observed 
that population evolved without dominance is larger than the other without 
dominance due the deaths in the former being bigger than the latter, since in 
these dominant positions the effective survival probability is the minimum value 
between the values present in both sets at that age. In the particular simulation 
shown, for ago = 2 it is noticed a decrease in the survival probability when the 
dominance is considered, since the ages 2, 13, 10, 15 were dominant positions. 

Figure 6 shows the time evolution of the total population of each age for 
sexual reproduction when the mutations are exclusively harmful. From this 
figure wo can notice that a stable population for ages a < 3 is obtained, in 
contrast to the original model in which there are no individuals alive older than 
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Figure 5: Survival probability for asexual reproduction (circles), sexual repro- 
duction with amin = IjCtmax = 17 and d = 1 (squares), 2 (triangles). Case (a) 
(top) and case (b) (bottom). 

age a > 2, even if beneficial mutation and also a deleterious mutation rate 5 
times smaller have been assumed. The result obtained here (Fig. 1), introducing 
sex in the original model, was found with the asexual Hotzcl model [15] only 
when mutations were allowed on a very small fraction F = 10% of the babies 
and also a minimum age of reproduction amin = 8 was considered. In our 
simulations, F = 100%, amin = 1 and amax = 17. 

4 Conclusions 

We have shown that main problem related to the Hotzcl model, which was its 
incapacity to treat populations with many age intervals, has been overcome 
by introducing recombination (with males and females) in this model, without 
any other assumptions. As well as, with the inclusion of sex in the model, the 
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Figure 6: The total population for ages 0, 1, 2, 3 (top to bottom, respectively) 
as a function of time for sexual reproduction when the mutations are exclusively 
harmful (ei = 0.0 and €2 = -0.20). 

population meltdown observed in the asexual version, when only deleterious 
mutations are considered, has been avoided. Moreover, in agreement with some 
earlier models, we have also obtained that the sexual reproduction (with males 
and females) produces better results than the asexual one. 
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